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The eigenfunctions and eigenenergies for a Dirac Hamiltonian with equal scalar and 
vector harmonic oscillator potentials are derived. Equal scalar and vector potentials 
may be applicable to the spectrum of an antinucleion imbedded in a nucleus. TYiaxial, 
axially deformed, and spherical oscillator potentials are considered. The spectrum 
has a spin symmetry for all cases and, for the spherical harmonic oscillator potential, 
a higher symmetry analogous to the SU(3) symmetry of the non-relativistic harmonic 
oscillator is discussed. 
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INTRODUCTION 


Recent theoretical investigations have suggested the possibility that the lifetime of an 


antinucleon inside of a nucleus could be significantly enhanced 


l|. The relativistic mean 


fields of antinucleons calculated in a self consistent Hartree approximation of a nuclear held 
theory 2, 3] indicate that the scalar Vs(f) and vector potentials Vy(r) are approximately 
equal. This implies that the antinucleon spectrum will have an approximate spin symmetry 
4, [a] • This conclusion is consistent with the fact that the nucleon spectrum has an approx¬ 
imate pseudospin symmetry 6, |7fl because the scalar and vector mean held potentials of a 
nucleon are approximately equal but opposite in sign and the vector potential changes sign 


under charge conjugation 


8]. In fact, the negative energy states of the nucleon do show 
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a strong spin symmetry [9|. Of course, whether such states could be observable can only 
be reliably estimated if the antinucleon annihilation potential is included in the mean held 
calculations. 

In this paper we shall solve for the eigenfunctions and eigenenergies of the triaxial, axially 
deformed, and spherical relativistic harmonic oscillator for equal scalar and vector poten¬ 
tials with the expectation that these results could be helpful in drawing conclusions about 
the feasibility of observing the spectrum of an antinucleon in a nuclear environment. The 


spherical relativistic harmonic oscillator with spin symmetry 


pseudospin symmetry 


lfl 


111 H 12, fl 14. [g 


and 


IT] have been studied previously, but in this paper we derive the 
eigenfunctions and eigenenergies for the triaxial and axially deformed harmonic oscillator as 
well. 


II. THE DIRAC HAMILTONIAN AND SPIN SYMMETRY 

The Dirac Hamiltonian, H , with an external scalar, Vs(r), and vector, Vv(r), potential 
is given by: 


H = a • p + P(M + V s (r)) + V v (r)) , 


( 1 ) 


where cc, /3 are the usual Dirac matrices, M is the nucleon mass and we set c — 1. The 
Dirac Hamiltonian is invariant under a SU( 2) algebra for two limits: Vs{f) = Vv(r) + C S and 
Vs(f) = —Vv(r) + C ps where C s ,C ps are constants 5]. The former limit has application to 


the spectrum of mesons for which the spin-orbit sp 


of an antinucleon in 
symmetry in nuclei 


itting is small 
. 


181 ] and for the spectrum 


he mean field of nucleons 8, 9]. The latter limit leads to pseudospin 
6|. This symmetry occurs independent of the shape of the nucleus: 


spherical, axial deformed, or triaxial. 
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A. Spin Symmetry Generators 


The generators for the spin 577(2) algebra, S q , which commute with the Dirac Hamilto¬ 
nian, [H s , S q ] =0, for the spin symmetry limit Vs(r) = Vv{r) + C s , are given by 5] 

\ 

( 2 ) 


S 9 = 


(l, (b 

( 

\ 0 ) 

V 


S q 0 


0 U n S n U n 


where s q = o q j2 are the usual spin generators, o q the Pauli matrices, and U p = ——— is the 


momentum-helicity unitary operator 


19|. Thus the operators 5* generate an SU(2) invariant 
symmetry of H s . Therefore each eigenstate of the Dirac Hamiltonian has a partner with the 
same energy, 

H„ <K Jr) = E,<K Jr) 


's J — ^k^k,/j,\ 


( 3 ) 


where k are the other quantum numbers and // = ±4 is the eigenvalue of S z , 


s z n,^ = ^HA?)- 

The eigenstates in the doublet will be connected by the generators 5±i, 

s± = 


( 4 ) 


'ji K 2 T ~^2 ± ^ ^ ±l( 


The fact that Dirac eigenfunctions belong to the spinor representation o 
given in Eqs. Q - leads to conditions on the Dirac amplitudes 20, 


the s 


21 


22 


( 5 ) 

pin SU(2), as 
, 23]. 


B. Dirac Eigenfunctions and Spin Symmetry 

The Dirac egenfunction can be written as a four dimensional vector 

/ + ^ \ 

9ljr) 


= 


9k,jr) 

iftiSf) 

\ kfkjr) ) 


( 6 ) 
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where g^Jj) are the “upper Dirac components” where + indicates spin up and — spin 


down and 44(4 are the “lower Dirac components” where + indicates spin up and — spin 


down. However spin symmetry imposes conditions on these eigenfunctions 
are derived from Equations 0 - © 


22 


2,3] which 


gti(r)=g i(f)=gk(r), 

*'2 ’ 2 

(4 

gt 44 = <444 = 0, 

K > 2 2 

( 8 ) 

4 + 44 = -4 44 = /44, 

^’2 ’ 2 

(9) 


(10) 


(11) 


Thus for spin symmetry the Dirac spin doublets are 


344 ^ 


( 0 ^ 

0 


#44 


> *t-*(4 = 


ifk(r) 


^4 + 44 

K, 2 

if ^/ 


v - ifk(f )) ) 


( 12 ) 


C. Second Order Differential Equation for the Eigenfunctions 

The Dirac Hamiltonian, ©, gives first order differential relations between gkif) and 
,4(4 and x (r). In the usual way we turn these equations into a second order differential 
relation for the upper component. In the limit of spin symmetry this second order equation 
becomes 

p 2 + 2 (. E k + M) V{r) - El + M 2 


9k(f) = 0 


(13) 
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where V s (f) = V(r) + V§, V v (r) = V(r) + V$, M = M + Vg, and E k = E k - V$. From Q 
and setting h= 1, the lower components become 


fk{r) = 


-1 


d 


(. M+E ) dz 


9k(r), 


f± (fA — 

k *2 ’ (M + E) K dx 


, d . d 

9k{r). 


(14) 

(15) 


These relations are consistent with the conditions on the eigenfunctions imposed by spin 
symmetry in Equations (Ell- 

Equation m is basically the energy-dependent Schrodinger equation without any spin 
dependence. Hence any potential, which can be solved analytically with the non-relativistic 
Schroedinger equation, is solvable in the spin limit, but the energy spectrum will be different 
because of the nonlinear dependence on the energy. Hence the Coulomb and harmonic 
oscillator are solvable. The Coulomb potential has been solved (for general scalar and 


vector potentials) 6] and applied to meson spectroscopy Il8| . More appropriate for nuclei 
is the harmonic oscillator. 


III. HARMONIC OSCILLATOR 

First we shall discuss the triaxial harmonic oscillator, then the axially symmetric har¬ 
monic oscillator, and finally the spherical harmonic oscillator. 

A. Triaxial Harmonic Oscillator 


For the potential V(r) = ^-(u)\ x 2 +u\ y 2 +lu 2 z 2 ), the second order differential equation 


m becomes 

3 2 d 2 d 2 


+ 7TW + 


(E ni ,n 2 ,n 3 + M) M(uj{ x 2 + u 2 y" + 1 X 3 z 2 ) + E 2 


dx 2 dy 2 dz 2 


n.i,n 2 ,n 3 


M 2 


9n\,n2,n 3 {+) 0- 

(16) 














6 


1. Eigenfunctions 

Introducing the product ansatz for the eigenfunction g ni ,n 2 ,n 3 (f) oc g ni ( x ) 9n 2 (y) 9n 3 {z), 
we derive the three equations 

d 2 


dxf 


— Xa + 2m + 1 


9n,(Xi) = 0, 


(17) 


where x\ = Ai x, x 2 = A 2 y, x 3 = X 3 z, and 


A,; = 


ni,n2,ns 


(^ni,n 2 ,n 3 + V) M Lijf 


M 2 = 2 £ (n,: + U 


(18) 

(19) 


The bound eigenstates are given by 


gni ,712,713 (f)=JV'(£? ni ,ri2,n3 ) ffm(a:i) SnaO^) ^ 3 (^ 3 ), 


where 


ShO^i) = 


A,; 


7T 2 ni np 


e 2 1 Hn.(xi), 


( 20 ) 


( 21 ) 


and H n .(xi) is the Hermite polynomial which means that g ni ,n 2 ,n 3 {f) has ni,n 2 ,n 3 nodes 
in the x, y, and z directions, respectively. M{E nim>n ^) is the normalization determined by 

JZc dx fZc d yfZc dz (\9nun2,n 3 (r)\ 2 + | fl, nm (f)\ 2 + | H = L 

From ED the lower components are 

fn\,n 2 ,n 3 (j) -Xf{E n i,n 2 ,n 3 ) 9n\(%l) 9n 2 {p^2) fn 3 +l (^ 3 ) 5 ( 22 ) 

= -^(^.^."s) (/m+lfcl) Snafo) T i g ni (x l) /na+l(® 2 )) f/n 3 (^ 3 ), (23) 


where 


frii+l (%i) 


A i 



H ni +i{xi) 

2 


- n H ni _ 1 (xi)'). 


(24) 
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Clearly the function f n+ \(x) is a polynomial of order n + 1. Evaluation for low n demon¬ 
strates that it has n + 1 nodes so we assume that it has n + 1 nodes for all n. This means 
that /ni,n 2 ,n 3 (^) has one more node in the ^-direction than g ni ,n 2 ,n 3 {f) and the same num¬ 
ber of nodes in the x and y directions as g ni ,n 2 ,n 3 {£) ■ On the other hand the amplitudes 
n:i ^(r) have the same number of nodes in the z-direction as g ni ,n 2 ,n 3 {^) ■ 

Using these amplitudes the normalization becomes 


N(E) 


2 (E + M) 


3 E + M 


(25) 


2. Eigenenergies 

From Equations (HIT A 11) the eigenvalue equation is 


\J £n\,n 2 ,n 3 T 1 (^-ni,n2,n 3 1 ) ^ ^ C> { n i T „)> 


(26) 


i =1 


where £ ni ,n 2 ,n 3 = Wl ^ 2 ’" 3 , Q = 2 and n* = 0,1,.... Thus both the eigenfunctions and 

the eigenenergies are independent of spin. 

This eigenvalue equation is solved on Mathematical 


Bni,n 2 ,n 3 d7 


1 4 

-B(Aii,n 2 ,n 3 ) + o + q d/4 \ 

O y ,712,713/ _ 


+ u°. 


(27) 


where 


and A 


n\ , 772,713 


-S( '^■711,712,713) 


— ELi ^ ( n * + 



A 2 

32 

Ani,n 2 ,n 3 T 1 

/ ^ 711,712,713 
/ 2 

27 





(28) 


The eigenvalues E nitn2>n3 are real for all values of ni,n 2 ,n 3 as long as Ci,Vg V are real. 
Although true it is not obvious because B(A nitn2iTl3 ) is not real for all A nin2 . n3 real. From 
d2H]), B(A ni) n 2 ,n 3 ) is clearly complex for A ni>n2tn3 < p. However, we now show analytically 


that E niri2ri3 will still be real even if B(A nit7l2jn3 ) is complex as long as \B(A nitTl2in 3 )| = 2/3. 















The imaginary part of E nitn2j1l3 is 

B(A ni ,n 2 ,n 3 ) + 1+9 B(An 1 ,n 2 ,n 3 ) - B (A nu n 2 ,n Z ) + I + 9 B(An 1 .n 2 ,n 3 ) 


Im E, 


ni,ri2,n3 


2 i 


Writing H(Au,n 2 ,n 3 ) = \B{A ni ^ n2in3 )\^ 

Im E nx ^ n2in3 f | B( y A ni n2 ri3 ) | 




9 |-B(^4n 1 ,n 2 ,n 3 ) | 


sin'll) , 


( 29 ) 


(30) 


and therefore E, 


is real if |-B(Aj 1 ,ra 2 ,ra 3 )| = 2/3 independent of ip. One can show 
numerically that |-B(A niin2in3 )| = 2/3 for all A ni n21l3 in the range from zero to -y/p, 0 < 
A ni ,n 2 ,n 3 < y/§- For A nijn 2 ! n 3 > §f, -B(An,n 2 ,n 3 ) is clearly real and hence E niiTl2>n 3 is real. 

The spectrum is non-linear in contrast to the non-relativistic harmonic oscillator. How¬ 
ever for small A r 


ni,n2,n3 


i ni,n2,ri3 


E ni ,n2,n 3 ~ M ( 1 -|- 3 +•'■) + Vy ( 31 ) 

and therefore the binding energy, E ni ^ n2 ^ n3 — M i=a J2i=i i n i + \) + Vy , in agreement with 
the non-relativistic harmonic oscillator. For large Aii,n 2 ,n 3 the spectrum goes as 

E ni ,n 2 ,n 3 ~ M (^4ni,n 2 ,n 3 + ~ + ' ' ' ) + Vy, (32) 

which, in lowest order, agrees with the spectrum for M 



B. Axially Symmetric Harmonic Oscillator 

For the axially symmetric harmonic oscillator uj\ — UJ 2 — cu_l and, hence, the potential 
depends only on and not the azimuthal angle (p , V(r) = ; 1 r l + z 2 ), where 

x = r _1 coscp, y = r± simp. This independence of the potential on p implies that the Dirac 
Hamiltionian is invariant under rotations about the £-axis, [H S 1 L Z ] = 0, where 
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and £ = r x p and and £ = U p £ U p . The Dirac eigenstates will then be an eigenfunction of 
L z and J z = L z + S z . 


Lz®N,n 3 ,A,»(r) — ^ 


h 4>K,n,,A,(.(0 = ^ A,aM- n = A + ft 


(34) 

(35) 


where N is the total harmonic oscillator quantum number and n 3 is the number of harmonic 
oscillator quanta in the z direction, and are discussed in more detail below. 


1. Eigenfunctions 

Since the potential has no 0 dependence the second order differential equation m sep 
arates into an equation for x 3 and an equation for p 

d 2 


dx 3 2 


- x-> ± 2n 3 + 1 


d 2 1 - 4 A 2 2 

[ap + ^~p p + 2n± + 2 


9n ,( 13 ) = 0, 

9n,(f>) = 0, 


(36) 

(37) 

(38) 


where p = \± r±, n± = 2 n p + A, and 


9N,n 3 ,A{ r ) 9n p {p ) 9n 3 {6 ^i 


A, = 


(E N , n3 + M) M cu]_ 


En,ti 3 ~ M 2 - 2 (A^_ ( n_ L + 1) + A3 ( n 3 + -)). 


(39) 

(40) 

(41) 


24], where 


The quantum numbers N, n±, n 3 , A are the “asymptotic” quantum numbers 
N = n± + n 3 is the total number of oscillator quanta, N = 0,1,..., n± is the number 
of oscillator quanta in the r^-direction, n± = 0,1,... IV, n 3 is the number of oscillator 
quanta in the ^-direction, n 3 = 0,1 ,N, and A is the angular momentum along the z-axis, 
A = ±n_ l, ±(t7j_ — 2) ■ ■ • ± 1, 0. The upper components of the eigenstates are given by 
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9N,n 3 ,k{j ) •A/"(-E l /V,n 3 ) A p 




n 


p- 


7 TA/ 7 T 2 " 3 (n p + A)! 713 ! 


o^+p 2 ) 

e ‘ 2 


/ <V) e‘ A * ff»,M,(42) 


where L^,\p 2 ) is the Laguerre polynomial which means that <Pv, B 3 ,a( 0 has n p ,n$ nodes 
in the re and z directions, respectively. A/"(i?v,n 3 ) is the normalization determined by 
fo* d( t> r r ± dr ± ^(^jv,n 3 ,A(r )| 2 + |/^ n 3 iAiP (r )| 2 + |/^ n 3 iA ^(f)| 2 ) = 1 and is the same 
function as given in H25D- These upper components do not depend on the orientation of the 
spin. 

The lower components are 


In. 


ri3,A 


r = 


A f(E. 


'N,ri 3 , A) 


n n 


M + En,u 3 


n (n p + A)! 


— £— A 

2 p A 


4 !V) 


fn 3 +i{x 3 ), (43) 


/ 


N, 723, A,— 



A/~(£^y in3 ) 2 / n p- A-i 

M + E N , n3 p V n(n p + A)\ 

(K + 1) A^+iV) + K + A ) A^“ 1} (p 2 )) c/ n3 (x 3 ), (44) 


fN,n 3 ,A,%(?) 


■Af(EN tri3 ) ^2 / 

M + E NjU3 p y tt (ra p + A)! 

+ 44 iV)) e 4 A + 1 ) < ^ > ^(* 3 ), 


P A+1 « +1) (P 2 ) 


(45) 


where g n 3 (x 3 ), fn 3 ( x 3 ) are defined in (I2T1) and (12-H) . Clearly the function /v,n 3 , a(A) has n 3 + l 
nodes in the ^-direction, one more node than the upper component, but the same number 
of nodes in the rj_ direction. The amplitude n A _i (r)) has the same number of nodes in 
the z - direction as the upper component. However, it has n p + 1 nodes in the rj_-direction 
for low n p and we assume it has n p + 1 node for all n p \ that is, one more node in the re¬ 
direction than the upper component. On the other hand the amplitude f~ n A i(r) has the 
same number of nodes in the ^-direction and the re-direction as the upper component. 
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2. Eigenenergies 

From Equations (HIT B 11) the eigenvalue equation is 

\J £jv,n 3 + 1 ( £n,u 3 — 1) — C± (n± + 1) + C 3 ('U 3 + —), (46) 

E 9 

where En, U3 = , Cj_ = — |^, and n p , n 3 — 0,1, .... Thus the eigenenergies not only have 

a degeneracy due to spin symmetry but they have an additional degeneracy in that they 
only depend on N and n .3 and not on A. 

The discussion about eigenergies is the same as for triaxial nuclei and the energy spectrum 
is given by 


En,u 3 = M 


B(A, 


n±,n 3 ) 


1 

3 


9 B(A n±iTl3 ) _ 


+ V°, 


where = C± (ni + 1) + C 3 (n 3 + |). 


(47) 


C. Spherical Harmonic Oscillator 


For a spherical harmonic oscillator uy = u and hence the potential depends only on the 
radial coordinate, r = \Jx 1 + y 2 + z 2 , and is independent of the polar angle, 9, z = r cos(9), 
as well as 0. The Dirac Hamiltonian will be invariant with respect to rotations about all 
three axes, [Lj, H s ] = 0 where 


L = 


(48) 


\ 0t J 

and hence invariant with respect to a SUl(2) x SU(2 ) group where SUl(2 ) is generated by 
the orbital angular momentum operators L. Since the total angular momentum, J = L + S, 
is also conserved, rather than using the four row basis for this eigenfunction, it is more 
convenient to introduce the spin function yy explicitly. The states that are a degenerate 
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doublet are then the states with j — t ± X and they have the two row form 22]: 




^ 9n r /(r) \YW(6,(j>) x\m ^ 


ifn r ,£,j ( r ) [^(M) X]£ 


(49) 


where ^ = £ ± 1 for j = £ ± Ym\d,(j)) is the spherical harmonic of order £, n r is the 

number of radial nodes of the upper amplitude, and [Y^\d, <f>) x]m i s the coupled amplitude 
^ C{Q■> 4>) Xfj,- Thus the spherical symmetry reduces the number of amplitudes in 
the doublet even further from four to three . 

The Dirac eigenstates will then be an eigenfunction of L ■ L, J ■ J, and J z , 


(r) = J (.j + !) KrAiM?)’ 

(50) 

■(r)=e(i + l) 

(51) 

n r ,l,j,M (r) = M ^ n r ,e,j,M (^)- 

(52) 


The differential equation for g nr /(r) becomes: 

(K - M 2 ) 


d 2 ^+1) , . 

dx 2 a; 2 


(E n + M) 


Mu} 2 


where E N = E N — Vy, M — M + VS, x = X r, and 


A = 


(. E n + M) Mu 2 


E% - M 2 = 2 A 2 (N + -) 


r g nr A r ) = 


(53) 


(54) 

(55) 


1. Eigenfunctions 


The solutions to this differential equation are well known and lead to the upper amplitudes 
of the eigenfunctions 

2A 3 n! £ ,(*+!) 


SW(r) = AT(E n/ ) 


r(£+n+ |) 


(56) 
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where J\f(E N ) is the normalization determined by j'^ d(p f* sin(6)d6 J 0 °° r 2 dr(\g nr ^(r)\ 2 
+ \fnr,t,j{r )| 2 ) — 1 an d is the same function as given in ( 1251 ) . Clearly g nr ,i( r ) has n r nodes 
in the radial direction. 

The lower components are 


f n r ,£,j =£— i (j ) 

A f{E N ) 


2 A 5 n\ 


M + E n V r<T + n r + §) 


"2 a* \(n r + 1) + (n r + £ + 


Ll h \x 2 )m 




2 A 5 n\ 


n r ,£,j=£+±V ) 


M(E N ) 

M + E]y V T(£ + Tl r + i) 


+A + AV)) 


(58) 


Clearly the function f nr ( j =i _ 1 (r) has n r +l nodes, one more node than the upper component. 
The amplitude f nr £ j =£+ i(r)) has the same number of nodes as the upper component. This 
agrees with the general theorem relating the number of radial nodes of the lower comonents 


to the number of radial nodes of the upper component 


25]. 


2. Energy Eigenvalues 


The eigenvalue equation is 


V£n + 1 {Sn - 1 ) = c (N + 1), ( 59 ) 

where £n = C = and N is the total oscillator quantum number, N = 2n + E = 
0,1,.... We note that there is not only a degeneracy due to spin symmetry but there 
is also the usual degeneracy of the non-relativistic harmonic oscillator; namely, that the 
energy depends only on the total harmonic oscillator quantum number and the states with 
orbital angular momentum £ = N, N — 2,..., 0 or 1 and angular momentum projection 
m = £, £ — are all degenerate. 

Again the eigenvalue is: 


En — M 


b[An) + 5 + 9 B(A n ) 


+ w 


(60) 
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and A n = C(N + |). 


D. Energy Spectrum 


In Figure Q] we plot the spherical harmonic oscillator Dirac binding enegies Ed = Ejy — M 
with Ejy given in Eq. EH, the solid curve, as a function of N. We chose the parameters to 
fit the lowest eigenenergies of the spectrum of an anti-proton outside of 16 0 in the relativistic 
mean field aprroximation 1] and they are C = 1.33, M = 252 MeV, and Vy = —677 MeV. 

The dashed curve is E D using the pertubation approximation of E^ given in Eq. (EH). 
The short-dashed curve is Ed using the asymptotic limit of E jv given in Eq. EH- Clearly 
the eigenenergies are in the relativistic asymptotic regime and not the linear regime of the 
non-relativistic harmonic oscillator. 

In Figure E| we plot the spherical harmonic oscillator excitation energies E^ = E^ — Eq 
for different N on the far left. Each level has a (N + 2) (.N + 1) degeneracy because of 
spin symmetry and because the allowed orbital angular momenta are £ = N, N — 2,... 0 
or 1 and the allowed orbital angular momentum projections are m = £, £ — 1, • • ■ — t. In 
the right of Figure El we plot the deformed excitation energies E* Nnz = E Nj713 — E 0 0 . The 
deformed excitation energies are staggered in groups for each N and each group contains 
the levels for 77.3 = 0,1,..., N with the excitation energy increasing with decreasing 77,3. 
The dimensionless oscillator strengths are determined by C 3 = C'j_ C 3 and assuming a 


deformation S = 0.33 which leads to C± = 1.49, C 3 = 1.05 


241. Each level has a 2 (IV—n 3 )-|-l 


degeneracy for (N — n 3 ) even and a 2 (IV — n 3 + 1) degeneracy for (.N — n 3 ) odd because 
of spin symmetry and because the allowed orbital angular momentum projections are A = 
±(IV — n 3 ), ±(N — n 3 — 2), • • • ± 1 or 0. The splitting of the levels within each N appears 
to be approximately linear with n 3 . 
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E. Relativistic Contribution 


The normalization A f(E) has the same functional form independent on whether the 
harmonic oscillator is triaxial, axially deformed, or spherical. This normalization has 
also been calculated independently by using f^dx f^dy f^dz (|/ / j h (F)| 2 + |/ / 7(b')| 2 ) = 
f^dx fZo d V f^dz g(r)* ^+^2 9 if) and we find agreement between the two different 
ways of calculating J\T(E). 

This also tells us that the probability of the lower component to the upper component is 
given by: 

f-o o dx IZ o d y IZ o dz (\fk,ZZ \ 2 + \fk,ZZ\ 2 ) _ E k -M 


Rl — 


(61) 


fZo dx fZLdyfZodz \9k(r)\ 2 2 (E k + M)' 

Thus for E k k, M the system is not very relativistic and the contribution of the lower 
components is small. For E k >> M, this ratio approaches |. For free particles this ratio 
approaches unity which indicates that the harmonic oscillator reduces the relativistic effect. 

In Figure 01 we plot this ratio for the spherical harmonic oscillator, Rjy, as a function of 
N. Even for the most bound states this probability is about 20 % and thus the antinucleon 
bound inside the nucleus is much more relativistic than a nucleon inside a nucleus for which 
this probability is about 1 %. 


IV. HIGHER ORDER SYMMETRY 


26]. This 


The non-relativistic spherical harmonic oscillator has an SU(3) symmetry 
symmetry is generated by the orbital angular momentum operators £ and the quadrupole 
operators 


a - lrrl (2) A 2 + 

( dm ~ I' 1 \m a NR ' \ 2 

A 


(62) 


NR 


where [rr]^ means coupled to angular momentum rank 2 and projection m and A nr — 
a [Wlo. This quadrupole operator is then a function of the dimensionless variable x^r = 
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Xnr f. These generators connect the degenerate states of the harmonic oscillator. 

The same degeneracy that appears in the non-relativistic spectrum appears in the rela¬ 
tivistic spectrum. The upper component of the relativistic eigenfunction given in Eq. m 
has the same form as the non-relativistic harmonic oscillator eigenfunction except that it 
is a function of the relativistic dimensionless variable x = Xr. Therefore if the generators 
are written in terms of the relativistic dimesionless variable they will connect the upper 
components of all the degenerate states of the relativistic harmonic oscillator in the same 
manner as the non-relavitistic quadrupole operator in Eq. © connects the degenerate 
eigenstates of the non-relativistic harmonic oscillator. Likewise, since the lower components 
are proportional to U p operating on the upper components, the dimensionless quadrupole 
operator transformed by U p connects the lower components of the degenerate states in the 
same manner as the non-relavitistic quadrupole operator in Eq. © connects the degener¬ 
ate eigenstates of the non-relativistic harmonic oscillator. However A depends on the energy 
( see Eq. (ED) and therefore the relativistic quadrupole generator is 


/ 


Qm 


\rr 


( 2 ) 


0 


\ 


^ 0 U p [rr]m } u p j 


lusJM(H s - V° + M) + \pp]W-JM(H a - Vy + M) (63) 


u 


since U p commutes with p. These generators along with the orbitial angular momentum, 
L , given by (TIKI) . connect the degenerate states with each other. Work on the algebra is in 
progress. For M —>• oo, Q m —> Q m , 

( r_i (2) „ \ 


Qm 


\rr\m 


0 


y o U p [rr]m u p J 
which forms an SU(3) algebra with L. 


u M + 


[PP\ 


( 2 ) 

m 


u M ’ 


(64) 


V. SUMMARY AND CONCLUSIONS 


We have derived the eigenfunctions and eigenenergies for a Dirac Hamiltonian with triax- 
ial, axially deformed, and spherical harmonic oscillator potentials and with equal scalar and 
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vector potentials. In all cases the Dirac Hamiltonian is invariant with respect to the 577(2) 
spin symmetry and thus the eigenenergies are independent of spin. For axially symmetric 
potentials the Dirac Hamiltonian is invariant with respect to a SU(2) x f/(l) group and 
the eigenenergies are degenerate with respect to the orbital angular momentum projection 
along the z-axis which generates the f/(l). For the spherical oscillator the eigenenergies 
are degenerate with respect to the orbital angular mometum and hence invariant under the 
SUl(2) x 577(2) group. These energies also have a higher degeneracy which is the same 
as the non-rclativistic harmonic oscillator; that is, they depend only on the total harmonic 
oscillator quantum number. The generators that connect these degenerate states have been 
derived but a larger symmetry group analogous to 577 (3) symmetry has yet to be identified. 
However, for infinite mass, the spherical relativistic harmonic oscillator is invariant with 
respect to an SU (3) x SU (2) group. 

The eigenenergies have the same functional form for triaxial, axially deformed, and spher¬ 
ical potentials and depend on one variable which is a linear combination of the oscillator 
quanta in a given direction weighted by the strength of the oscillator potential in that 
direction. The spectrum is infinite and the eigenenergies and are linear in the harmonic os¬ 
cillator quanta for small oscillator strength but increase slower than linear for large oscillator 
strength. 
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FIG. 1: The Dirac binding energies, Ejj, for the spherical harmonic oscillator as a function of 
N. The exact enegies are the solid line, the perturbation appoximation is the dashed line, and the 
asymptotic approximation is the short dashed line. 
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FIG. 2: On the left the excitation energies, E* N , for the spherical harmonic oscillator are plotted 
for each N. Each level has a (N + 2) (N + 1) degeneracy because of spin symmetry and because 
the allowed orbital angular momenta are l = IV, IV — 2,... 0 or 1 and the allowd orbital angular 
momentum projections are m = £, £ — 1, ■ ■ ■ — i. On the right the excitation energies, E*. for the 
deformed harmonic oscillator are plotted in staggered groups for each N. Each group contains the 
levels for = 0,1 ,,N. Each level has a2 (N — 713 ) + 1 degeneracy for (N — n 3 ) even and a 
2 (IV — ri 3 + 1) degeneracy for (N — n^) odd because of spin symmetry and because the allowed 
orbital angular momentum projections are A = ±(IV — 713 ), it (IV — 77,3 — 2), • • • it 1 or 0. 
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FIG. 3: The ratio of the probability of the lower components to the upper components, Rn , for 
the spherical harmonic oscillator as a function of N. 



